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M-P 逆 在 线性 随机 方程 中 的 应 用 ” 


姚 落 根 22， 欧 ELS MARRE 
(1- 湖南 师范 大 学 数学 与 计算 机 科学 学 院 ， 长 沙 410081; 2- 湖南 商学 院 信息 学 院 ， 长 沙 410205) 
摘 XE 本 文 证 明了 可 料 过 程 的 M-P 逆 仍 然 可 料 。 利 用 M-P 逆 讨 论 了 线性 随机 方程 可 料 解 的 结构 ， 由 此 
求 出 了 扩散 模型 中 的 所 有 等 价款 测度 ， 并 给 出 了 最 小 逆 炽 测度 。 
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1 引言 


在 金融 数学 及 其 他 领域 研究 中 ， 经 常 遇 到 如 下 线性 随机 方程 
oOY(G = C(D，t> 0， (1) 


其 中 5，C 分 别 是 给 定 概率 空间 (N, F, Fe P) ES R^" (EM R EENE. RAT 
要 得 到 (1) 的 所 有 可 料 解 。(1) 有 解 (未 必 可 料 ) 的 充 要 条 件 是 


ranke(t) = rank(c(t), C(t)), vt20 


几乎 处 处 成 立 。 本 文 利 用 高 等 代数 中 的 M-P 道 来 求 这 个 方程 (rc 未 必 几 乎 处 处 满 秩 ) 的 可 料 解 。 
定义 1 REAR”, WREX e R 满足 如 下 四 个 方程 


AXA= A, XAX-X, (AXY- AX, (XA) = XA, 


则 称 半 为 4 的 M-P 逆 ， 记 为 4+。 

Moore 和 Penrose 已 证明， 对 任意 实 和 矩阵 A， 其 M-P 道 是 唯一 存在 的 。 下 面 以 引 理 的 形式 
给 出 将 要 用 到 的 性 质 。 

9381 设 Ae Rdx*， 则 

1) (A)* - (A*y, A* = (A A)* A = A(AA)+, (AA)+ = A*(A*Y's 

2) ”如 果 线 性 方程 组 Ar =b AN, WX TBBUSERDIAWT-—A*tb4(—A*A), 
其 中 z 是 任意 的 4d 维 列 向 量 。 

对 R4 值 ， 可 料 随机 过 程 M(t) = (Mi(t).--- ,M2(t))， 若 M(t) e M2,,.( 初 值 为 0 的 局 部 平 
方 可 积 蒜 )， 则 其 二 次 变 差 过 程 为 


(Mi(t), M(H) =- (M(t), Malt) 


(M(t) = 
(Malt), Mi) -> (Malt), Malt)) 
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X [1] WEHA T (M(t)) 的 M-P 逆 仍然 可 料 ， 并 且 以 一 个 反例 来 说 明 ，(M 人 Gy 的 其 他 广义 逆 不 一 定 
保持 可 料 性 。 受 此 启发 ， 本 文 证 明了 任意 取 值 于 Rdx"， 可 料 随机 过 程 的 M-P 逆 仍 然 可 料 ， 并 
由 此 求 得 随机 方程 (1) 的 所 有 可 料 解 。 

2 R”*" 值 可 料 随机 过 程 的 M-P 3 


本 节 考 虑 随机 过 程 oc 的 M-P 逆 。 令 7 = ranko(t)， 则 7 可 料 趾 。 按 满 秩 分 解 定理 ， 对 (wt) € 
(r =k); FEITE Felt) € RIF, Gilt) e Rx*xn， 使 得 


a(t) = Ip F(t) Ge(t), 


其 中 rankF(t) = rankGx(t) = ko FTE 


d^n 


c(t) = >》 Ip FK(0)G&(t). 
k=0 
定义 g(t)+ 如 下 
d^n 
o(t)* = V Irse) Gr lt) [GG (0/] ^ [FO Fkt) Filt). (2) 
k=0 


定理 1 Rot) TE, BDfz(2)5EXH)o(t)* J&o(t) H] M-PXÉ, H.o(t)* aE 

证 明 容易 验证 o(t)+ 满足 定义 1 的 四 个 方程 。 由 M-P XKE, n Tf o(2)* XE olt) 的 M- 
P 逆 。 还 需 证 明 o(t)+ 可 料 。 因 为 有 (如 ，Gx(t) 是 通过 可 料 过 程 o(t) 的 连续 变换 得 到 的 ， 
MF), GE) TEH FÆ, [Gk (0Gi(£)] t F [Fe (0) F(E] I EE Bir, 的 可 料 性 ， 即 
可 推出 c( 罗 + 可 料 。 

接 下 来 考虑 随机 方程 (1) 的 解 (假定 其 存在 解 )。 首 先 需 区 分 可 料 解 与 非 可 料 解 。 

定义 2 FRAMI O = {6B(t);t > 0} 为 方程 (1) 的 一 般 解 ， 如 果 


P(fc(t)O(t) = C(t), Yt > 0})=1. (3) 


如 果 日 还 是 可 料 的 ， 则 称 为 严格 解 。 

显然 ， 严 格 解 一 定 是 一 般 解 ， 但 一 般 解 不 一 定 是 严格 解 。 

定理 2 ”随机 方程 (1) 存在 严格 解 的 充分 必要 条 件 是 它 存在 一 般 解 。 

证 明 如 果 Y(t) 是 随机 方程 (1) 的 严格 解 ，Y(t) 必 是 方程 (1) 的 一 般 解 。 现 在 设 Y(t) 是 随 
机 方程 (1) 的 任意 一 个 一 般 解 ， 令 O(t) = o(t)+C(t)。 则 由 定理 1，Bo(t) 是 可 料 过 程 。 再 由 


c(t)Oo(t) = o(t)o(t)* C(t) = e(t)o(t)* o(t)Y (t) = e(t)Y(t) = C(t). 


故 Bo(t) 是 随机 方程 (1) 的 严格 解 。 
下 面 讨论 随机 方程 (1) 解 的 结构 。 令 


Alo) = (Y (t) : Y (t) 是 取 值 于 R” 的 可 料 过 程 , Ho(t)Y (t) 20, t 2 0, a.s.P-), 
Blo) = {(In 一 o(t)to(t))2Z(t) : Z(t) ZERAT R” 的 可 料 过 程 .}， 
其 中 ,是 nn 阶 单位 阵 。 因 为 Y(t) 三 0e A(o)， 故 A(o) 不 是 空 集 。 
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定理 3 A(o) = B(o). 

证 明 设 Y(t) e A(o). WE (I, — o(t)*o(t))Y (t) = Y(t), 知 A(o) C Blo). gi Z(t) Æ 
取 值 于 R” 的 任意 可 料 过 程 。 因 为 

o(t)(In — o(t)to(t))2Z(t) = e(0Z(t) - a(t)a(t}ta(t)2Z(t) = e(t)Z(t) — a(t)Z(t) = 0, 

故 有 A(o) 2 Blo). EEE. 

定理 4 O(t) (1) 的 严格 解 的 充 要 条 件 是 存在 Y(t) € A(o)， 使 得 9@(t) = (t) +Y (t). 

证 明 i e(t 是 随机 方程 (1) 的 严格 解 ， 则 e(t) 可 料 。 令 Y(t) = eq) 一 GQo(t)， 显 
R YATE. HE 

a(t)Y (t) = ao(t)O(t) ~ e(t)Oo(t) = 0. 


MY (t) e A(o). 反之, Y(t) e Alo) MOG = 6o(b --Y(t) TH H 
ce = e()Y (t) + o(t)Oo(t) = o(t)Go(t) = C(t). 
H O(t) 是 随机 方程 (1) 的 严格 解 。 
由 定理 3 和 定理 4， 随 机 方程 (1) 的 严格 解 的 通 解 具有 如 下 形式 
e(t) = Oolt) + (In — o(t)* e(t) Z(t), (4) 


其 中 Z() 是 取 值 于 Ra 的 任意 可 料 过 程 。 


3 M-P 逆 在 扩散 模型 中 的 应 用 
考虑 如 下 扩散 模型 ， 市 场 中 有 4 十 1 种 资产 ， 其 价格 过 程 满足 如 下 微分 方程 


| dSo(t) =r#S0(t)dt, So(0)=1, O«t«T, 
(5) 


d5;(t) = ai(t)Si(t)dt + 5i(t) » cij(t)dW;(t), Si(0)= Si, i=1,...,d, 
je 


其 中 
W = (W(t) = (Wi Walt), ---, Ws). Fu O<t<T] 


是 定义 在 某 带 ce- 代数 流 的 概率 空间 (Q, F, Fa, P) ERU b ER BH3S JJ. {Fi} 是 由 自然 -代数 
流 (EN) 在 概率 已 下 的 扩张 形成 的 。 假 定 平均 收益 率 过 程 (t) = (a(t), a(t), --- ,aa(0) F 
唯一 强 解 所 需要 的 条 件 (如 Lipschitz 条 件 }。 为 简单 起 见 ， 记 上 面 的 市 场 模型 为 WM(r,a,c)。 
WCU) = a(t) -rt 了 ， 这 里 了 是 每 个 分 量 都 为 1 的 d 维 矩阵 。Q 是 等 价 蒜 测度 的 充 要 条 件 
是 存在 可 料 过 程 @(t),，P( 广 9(w)|j?du < +o0) = 1， 使 得 
dQ 


dP 


Se «( = I O(u)dW(u)), (6) 
HE Ol) = (6B1(t),… ,Bn(t)) 满足 如 下 的 随机 方程 


c(t)X(t) -C(t, P-as. te(0,T]. (7) 
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定义 3 称 取 值 于 及 "的 随机 过 程 9 = {6(t), 0 <t < TAARA, WR OTH, 
且 满 足 (7)。 

为 了 求 出 一 个 具体 的 风险 溢价 过 程 ， 通 常 对 c(t) 作出 许多 限制 ， 目 的 都 是 保证 [c(t)o (t)! 
RÉ [o(t)a(0]-? FE. FIECBUTERIRI M-P 逆 来 解决 这 个 问题 。 如 果 记 也 为 M(rua,o) 中 等 价 
BWES, $ 

T 
A= [zo : Z(t) © d ERES, HP (f ll — e(u)*oe(u)) Z(u)|?du < +oo) - 2: 
0 
定理 5 如果 Bo(t) 满 足 


T 
BC lOo(u)l?du < oo) = 1. 
则 Q e P 的 充 要 条 件 是 存在 Z EA, TE 
= |a -«(- J í Bo(waw (u))e( - f ERE e(u)*a(u)) Z(u)dW (u)). 
证 明 ”首先 注意 到 风险 溢价 过 程 6(t) 具有 形式 (4)， 且 有 
(8o(t))' Us — e(t)*o(0)Z(t) = (HTC n — o(t)*o(0)Z(t) 
= (c(t)*e(t)Go(0) (In — c(t)*o(0))Z(t) = 0. 
所 以 有 
le = [Go(t) + Qs — e(t)*e()Z(Ol'Ieo(t) + Un — o(t)*o(0)Z(0] 
= [OH + IG — e(t) *e()Z CI. 
在 定理 的 条 件 下 ， 显 然 
T T 
P( T lO(u)l?du < +0) =1 e P( f ln — olu)to(u)Z(u)lPdu < 49) = 1. 


即 Z e A。 再 由 (6)， 
E E = exp { — [ O(u)dW (u) 一 J ledolsav} 


= «(- I l Oo(u)dW (u))e( — f Te e(u)*e(u)) Z(u)aW (u)). 
idis Z(t) 三 0 对 应 的 拷 测 度 为 Qo， 即 


d t 
ae |^-«(- f ev(u)dW(u)). 
由 定理 5， 蒜 测度 之 间 的 差别 ， 仅 仅 在 于 它们 对 应 的 可 料 过 程 Z 不 同 。 于 是 ， 选 择 哪 一 个 著 测 
度 作为 定价 测度 的 问题 就 转化 为 求 其 对 应 的 Z。 这 样 一 来 ， 问 题 就 简化 了 许多 。 下 面 的 例子 说 
明了 这 个 问题 。 
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例 1 在 模型 M(ra,o)'P. RIMAR (the minimal reverse entropy martingale mea- 
sure) JJ Qo. 
T idt 8 du JE Ze ee C A ja] A 


T 
E»(- m3) = ie/ testo 1a. ol) tal) zlau) 


IV 


T 
5E (f. leotada). 
即 Qo 是 最 小 着 病 半 测度 。 
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Application of Moore-Penrose Inverse in Linear Stochastic Equation 
YAO Luo-gen"?, OU Huit, YANG Xiang-qun! 


(1- College of Mathematics and Computer Science, Hunan Normal University, Changsha 410081; 
2- College of Information, Hunan University of Commerce, Changsha 410205) 


Abstract: This paper proves that the M-P inverse predictable process is still predictable. By employing 
the M-P inverse, we discussed the structure of predictable solutions of the linear Stochastic equation. 
All equivalent martingale measures and the minimal reverse entropy martingale measure in a diffusion 
model are obtained. 
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